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We study theoretically the collective orbital excitation 
named orbital wave in the orbital ordered manganites. The 
dispersion relation of the orbital wave is affected by the static 
spin structure through the coupling between spin and orbital 
degrees of freedom. As a probe to detect the dispersion rela- 
tion, we propose two possible methods by utilizing resonant 
inelastic x-ray scattering. The transition probability of the or- 
bital wave scattering is formulated, and the momentum and 
polarization dependences of the structure factor are calculated 
in several types of the orbital and spin structures. The elas- 
tic x-ray scattering in the L-edge case to observe the orbital 
ordering is also discussed. 

71.10.-w, 71.90.+q, 75.90.+W, 78.70.Ck 



I. INTRODUCTION 

It is widely accepted that one of the important in-, 
gredients for the colossal magnetoresistance (CMR)trcl 
and various complex phenomena observed in perovskite 
manganites is the orbital degree of freedom in Mn 3+ ion. 
Although the long range orbital ordering and its impli- 
cations oLthe magnetic interaction were discussed long 
time agoJjLj it has been recognized as a hidden degree 
of freedom, because the observation technique has been 
limited.EI Recently, the anomalous x-ray scattering (the 
resonant elastic x-ray scattering) shed light on the sub- 
ject, that is, by using the method the orbital ordering 
in Lao.5Sri.5Mn04 was directly observed.EI The alternate 
orbital alignment was confirmed by observation of the 
forbidden reflection originated from orbital dependence 
of the atomic scattering factor. In a moment, this ex- 
perimental technique was recognized as a powerful probe 
to detect the orbital structure-and applied to several un- 
doped and doped manganites. E3H3 Though the intensive 
experimental and theoretical studies ,Ej~t3 it certainly un- 
covers roles of the orbital ordering on the electronic and 
lattice structures in manganites. 

On the other hand, the dynamics of the orbital degree 
of freedom is still far from our understanding . In orbital 
ordered insulators, the collective excitation for the orbital 
degrees was first investigated by Cyrot et al. in the model 
where the orbital space was assumed to be isotropicLa 
It was termed orbital wave in analogy with spin wave 
in a magnetic ordered state. For actual compounds, 
the present authors calculated the dispersion relation of 



the orbital wave in the (Sc^a _ r 2/3d3j,2_ r 2)-type orbital 
ordered state with A(layer)-type antiferromganetic spin 
structure observed in LaMn03.E3 When such a collective 
excitation in the orbital degree exists, it probably affects 
the thermodynamics and the dispersion relations of spin 
waveliS and phonon. Furthermore, the orbital excitation 
in the metallic phase seems to play .an essential role in 
the transport and optical propertiest3 and CMR. As a 
probe to observe the orbital excitation, Inane et al. the- 
oretically examined the Raman scatteringpj where the 
excitation process is analogous of that in the magnon 
Raman scattering. In this method, the information of 
the T-point and the density of state of the orbital wave 
are obtained. 

In this paper, we study theoretically the orbital wave in 
orbital ordered manganites and propose a method to ob- 
serve it by the resonant inelastic x-ray scattering (RIXS). 
RIXS is rapidly developed through the recent progress of 
the synchrotron radiation sorcecil and is_applied to sev- 
eral highly correlated electron systems.Eju 3 ] One of the 
advantages of this method on the subject is that the wave 
length of the x-ray is comparable to the lattice constant, 
unlike the conventional light. Therefore, the dispersion 
relation of the orbital wave in the wide range of the Bril- 
loun zone is detectable by the method. The following 
two different processes in RIXS are proposed: the inci- 
dent photon energy is tuned at 1) Mn 3+ L-edge and 2) 
Mn 3+ K-edge. We formulate the transition probability 
in both types of the scattering and the structure factor is 
calculated as a function of momentum and polarization 
in several orbital orderings. 

In Sec. II, the dispersion relation of the orbital wave is 
calculated. We emphasize the strong correlation between 
the static spin structure and the orbital wave. In Sec. 
Ill, a mechanism to observe the orbital wave by RIXS 
is introduced. The general formulas of the scattering 
intensity and the numerical results of the structure factor 
are presented. Sec. V is devoted to the discussion and 
summary. 



II. ORBITAL WAVE 

In order to calculate the collective excitations in an or- 
bital ordered insulator, we start with the Hamiltonian de- 
scribing the low energy electronic structure in perovskite 
manganites. The cubic lattice consisting of Mn ions is 
considered and two e g orbitals and a ti g localized spin 
are introduced in each site. The Coulomb interactions 
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between e g electrons (U, U' , J) and the Hund coupling 
(Jh) between e g and t g spins are considered. Because 
the Coulomb interaction gives the largest energy among 
the relevant parameters, the Hamiltonian is derived by 
excluding the doubly occupied e g states as follows, 

H 3 d = H e -e + H e . t + H t - t . (1) 

The detailed derivation is presented in Ref. [l7| The first 
term describes the interaction between nearest neighbor- 
ing spins and orbitals in e g orbitals, 

H e . e = - 2J a + Si • 4) (\ - 

(ij) 

- 2J i E(i - §i ■ § i) (I + + Ai d ■ ^ 

where a relation U = U' + J is assumed with U, U' , and 
J being the intra- and inter-orbital Coulomb interactions 
and the exchange interaction in e g orbitals, respectively. 
Ji = tl/(U'-J) and J 2 = tl/(U'+J+2J H ) with J x > J 2 . 

=[Ti Z ,Ti x ] is a doublet of the orbital pseudo-spin 

operator: Tj = (1/2) X) CT77 ' A-yA^i-f^ii'" and S ije is 
the spin operator with S — 1/2 for an e g electron of 
orbital jg. fy and A^ are defined by 

l(l + cosni?f sinn;^ \ .. 
z y smri;— 1 — cosn;— y 

and 

A i:j = ^tMijiij + H.c. , (4) 

respectively, where Mu+i = (cosn/^, — sinnj^)* with 
(n x ,n y ,n z ) — (1,2,3). The forms of and My in the 
above are obtained by the Slater-Koster formulas for the 
transfer intensity tj^ between site i with orbital 7 and 
site j with 7'. It is worth noting that the orbital part of 
H e - e is represented by T z and T x , since t/J is rewritten 
by the rotation matrix of the (T z , T^-plane and the pro- 
jection operator for 3d3 Z 2_ r 2 orbital. The first term in 
Eq.(H) favors the ferromagnetic spin structure with the 
antifcrro-type orbital ordering where two kinds of orbital 
sublattice exist. On the other hand, the second one fa- 
vors the antiferromagnetic structure with the ferro-type 
orbital ordering where occupied orbital is uniform. The 
sum of the second and third terms in Eq. (^) is given by 

H e -t + Ht-t = — Jh St 2g i • Si le 

+ Jaf ^2 St w ' St 2g j , (5) 
(«> 

where Jaf is the antiferromagnetic superexchange inter- 
action between nearest neighboring t 2g spins St 2 % with 
5 = 3/2. 



We first study the spin and orbital ordered structures 
at zero temperature in the mean field approximation. It 
is assumed that in each Mn site, one of the e g orbitals is 
occupied by an electron. For both spin and orbital struc- 
tures, we introduce four-types of the ordering: ferro-type 
(F-type) where spins or orbitals in all sites are parallel, 
and three kinds of antiferro-type, that is, layer-type (A- 
type), rod-type (C-type), and NaCl-type (G-type) where 
two kinds of spin or orbital sublattice exist. The rotating 
frame is introduced in the orbital space and the orbital 
state is described by an angle 9 in the (T z , T x ) -plane. 
In the scheme, the occupied orbital is represented by 
|3d 7e ) = cos(9/2)\3d 3z 2_ r 2) — sm(9/2)\3d x 2_ y 2). As the 
order parameters, (S z ), (St 2gZ ) an d (T(9) z ) — cos9(T z } — 
sm9(T x ) with (St 2 z ) — 3{S Z ) are introduced. 

In Fig. 1, the mean field phase diagram is presented. 
A sequential change in the spin structure, that is, F — > 
A — > C — > G with increasing J 2 and Jaf is caused by 
an enhancement of the antiferromagnetic interaction be- 
tween spins. As for the orbital structure in the spin-F and 
spin-G phases, the orbital-G with (9a/9b) = {9a/Sa + ^) 
for any 9 A and the orbiteLCwith {9 A /9 B ) = i^/^f) are 
the mean field solutionsQoE£l with 9a and 9b being the 
angles in A and B orbital sublattices, respectively. In 
the orbital-G case, especially, the orbital space becomes 
isotropic in the (T z , X^-plane, although the orbital part 
in the Hamiltonian Eq. (|^) is not written by a simple vec- 
tor product. This result is originated from 1) the relation, 

E *£& = !*o&n< , (6) 

l=x,y lZ 

due to the cubic symmetry of the lattice and 2) the cu- 
bic symmetry of the spin structure, that is, spin-F and 
-G. Because of these two conditions, the orbital part in 
H e - e becomes diagonal and Aij term vanishes. On the 
other hand, in spin-A and -C phases, the orbital state 
is uniquely determined. In spin-A, the orbital structure 
in the mean filed theory is orbital-C and -G types with 
[9 a/ 9 b) — {9 a I ~ 9 a) where 9 a is determined by the 
equation 9a — cos _1 (2J 2 /(5Ji — J 2 )). It includes the 
(7y=(3(i3 Z 2„ r 2 — Zd x 2_ y 2) / ^(2,d 3z 2_ r 2 +id x 2_ y 2)) orbital 

ordep-,at| J 2 = and (3d z 2_ a ,2/3(ij,2_ z 2) order at J 2 — 
Ji.tdoEj The (6a/ ~ 9 A )-type orbial order in spin-A is 
consistent to the experimental results of the polarization, 
dependence of the resonant elastic x-ray scatteringjlj'til 
although the value of 9a seems to be modified by the 
coupling with lattice. In spin-C, the mean field solu- 
tion is orbital-G type with (0/n) corresponding to the 
(3c?3 Z 2_ r 2/3(i 2 ;2_j ( 2). In both cases, the spin structure 
breaks the cubic symmetry in the system and causes the 
anisotropic interaction between nearest neighboring or- 
bital pseudo-spins. This noticeable feature in the inter- 
play between spin and orbital degrees clearly reflects on 
the dispersion relation of the orbital excitation, as will 
be discussed below. 

We next study the orbital excitation in the orbital or- 
dered states obtained in the mean field theory. The col- 
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lective excitations in the spin and orbital degrees of free- 
dom are obtained by utiligkig-,the conventional Holstein- 
Primakoff approximationfiS ea The dispersion relation in 
the spin-F case is shown in Fig. 2 where the types of the 
orbital order are chosen as 9 a — Trn/6 = 9 B — 7r. The 
face-centered cubic (fee) lattice is adopted as a unit cell 
which includes two Mn ions. The analytic form of the 
dispersion relation and the eigen-operator are given by 



with 



a 



4 ±} = 3^/a(a±P) , 



l—x,y,z 



l 



/3 = 2Ji^- \ cos(aki) (t u+ i) 22 , 

L—x,y,z 



(7) 



(8) 



(9) 



and 



cosh^ ±) -^=(a fc ± b k ) 

- Sinh ^ ±) 7^( a -fc ±6 U): 



(10) 



with 29^ = tanh _1 (±/3/(2a±/3)), respectively, a& and 
bk are the Fourier transform of the Holstein-Primakoff 
bosons for the two orbital sublattices defined in the ro- 
tating frame. The (+)-mode is the gapless mode and its 
eigen-operator includes the components Tax — T Bx and 
Tatj + T By . At the T-point, the cc-component of the stag- 
gered orbital order parameter ({J2ieA ^fe - J2teB T ix) 2 ) 
diverges since the orbital space is isotropic. With in- 
creasing 9a, the stiffness of the orbital wave along the 
r — X direction becomes weak and at 9a — 47r/6 
(3d,3 X 2_ r 2 /3d y 2_ z 2) it becomes fiat. It is originated from 
{tu+x)22 — in Eq.(Q) since the electron transfer between 
id y 2_ z 2 orbitals along the x direction vanishes. The sim- 
ilar gapless mode and the changes of the stiffness with 
changing types of the orbitals are also obtained in the 
spin-G case. The isotropic nature in the orbital space 
in spin-F case is seen by the Fourier transforming of the 
Hamltonian (Eq.(||)) as 

H e - e = -2Ji^(J - mMmik)) , (ii) 

k l—z,x 

in the spin-F case@ Ti(k) is the orbital pseudo-spin op- 
erator which diagonalizes the matrix Tjy and 



Ti(k) = 2((c x +c y + c z ) 



±(<2 



CxCy CyC z 



Cz Cx ) 



1/2^ 



(12) 



where + and — are for I — z and x, respectively, and 
q = cos(afc;). At the T- and i?-points, r x = t z , so that 
the orbital system becomes isotropic. 



The orbital excitation in the spin-A case is qualita- 
tively different from that in the spin-F case as shown 
in Fig. 3. The unit cell with four Mn ions is adopted, 
so that there exist four modes. Two of them have the 
dispersion relations, 



,(+±) 



- (2a xy + a zz ) [2{a xy + fl xy ) + (a zz ± (3 ZZ )] 



(13) 



and the other two cu^ ^ are given by replacing (a xy 

ftxy 

) in Eq.Jl3|) by ( 

&xy 

- (3 xy ). Here, 



ai,; 
and 



— 2Ji + Ki m j (T; m )ii 

+ 2J 2 Q -K lm y{T lm )n+2(M lm ) 1 ) , (14) 



2 Ji ( - + Ki , 



2J, 



(i 



(7, 



lm 22 



(15) 



for (lm) = (xy) or (zz). In these formu- 
las, Ai m = (l/2)E n= i, m 4+n for A = t and 



7/r, 



(1/2) E, 



=1 m r «+« cos(afc„) and K h 



is the spin correlation function given by Ki m = 
C 1 / 2 ) Y, n =i,m(( S i)z( S i+")z)- Among the four modes, 
(++)-mode is the lowest one and its eigen-operator is 
expressed as 

<t> ( k ++) = cosh6^ ++) i(aife + a 2k + h k + 6 2fc ) 

+ S mb4 ++) i(al_ fc + 4_ k + b\_ k + b\_ k ) , (16) 

where the subscripts 1 and 2 imply the spin sublattices. 
The operator includes the components Ta iX + Ta 2 x — 
T Blx - T B2X and T Aiy + T Am + T Biy + T B2y and it corre- 
sponds to (+)-modc in the spin-F case. At the T-point, 
the energy is obtained as 



,(++) - 



J k=r 



= - V (2a xy 



,)A {++) 



(17) 



with 



A (++) 



(lm) = {xy),(zz) 



Zlr 



2J 1 (-+K lm } +2J 2 (- 

( T im)22) 

2(M lm )i 



Ku 



l 



2J 2 [--K l 



(18) 



where z xy = 2 and z zz = 1. This mode is gapful in 
contrast with the spin-F case. The right hand side in 
Eq. ( |l8| ) is represented by a product of the spin cor- 
relation function and the anisotropic transfer intensity. 
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Therefore, the origin of the gap is attributed to the fact: 
1) The orbital space is anisotropic, i.e., (r)n 7^ ( T )22 
and Aij ^ 0, which originate from the anisotropic trans- 
fer intensity tj^ ^ S 71 >. 2) The spin correlation func- 
tions in the xy plane and in the z direction are different 
{K xy ^ K zz ). The gapful excitation is also obtained in 
the spin-C case. It is concluded that the anisotropic spin 
structure breaks the cubic symmetry in the system and 
causes the anisotropy in the orbital space, which causes 
the gap in the orbital wave. 



determines the direction between i and j sites, unlike in 
the present case. In the case of manganites, the charge 
transfer from the nearest neighboring Mn sites occurs 
more easily than cuprates and nickelates, since the intra- 
site Coulomb interaction (U) is smaller and the energy 
gain of the Hund coupling in the intermediate states (3d 5 
states) is larger in the present case£il In the L-edge case 
the one-orbital wave excitation is only brought about. 
On the other hand, in the K-edge case, both one- and 
two-orbital wave excitations occur. 



III. RESONANT INELASTIC X-RAY 
SCATTERING 

In this section, we propose a method to detect the 
orbital wave by RIXS. The following two excitation pro- 
cesses are considered: 1) the incident energy is tuned 
at Mn 3+ L-edge and 2) at Mn 3+ K-edge. Hereafter, the 
former and latter processes are termed L-edge and K- 
edge cases, respectively. The schematic processes are pre- 
sented in Fig. 4. In the L-edge case, the incident x-ray 
excites an electron from Mn 2p orbital to the unoccupied 
Mn 3d orbital, and then one of the 3d electrons fills the 
core hole with emitting the x-ray. When the occupied 
orbitals are different between the initial and final states, 
the orbital excitation is brought about and the excitation 
propagates through the interaction between the nearest 
neighboring pseudo-spins. This process is denoted by, 



\3di) + hv^ \3d\3dl,2p) -> |3dL) + hv' 



(19) 



where 2p implies that one hole occupies the Mn 2p or- 
bital. In the K-edge case, the excitation occurs through 
the screening process. At first, the incident x-ray excites 
an electron from Mn Is orbital to the unoccupied Mn Ap 
one at i site. At the time, an O 2p electron comes from 
one of the nearest neighbor oxygen sites to the Mn site 
and screens the potential due to the Mn Is core hole. 
This screening process is confirmed by the theoretical 
calculation in the small cluster rSiid the x-ray absorption 
spectroscopy experiments. EjEjo Due to the hybridiza- 
tion between the O 2p and Mn 3d orbitals, this state 
strongly mixes with the state where an e g orbital in one 
of the nearest neighbor Mn sites (j site) is empty. When 
the electron in the 4p orbital is relaxed to the Is orbital 
with emitting x-ray, one of the 3d electrons in i site comes 
back to j site and the orbital excitation is brought about. 
This is represented by 



|3dJ 7i 34 



31 3 I 



hv - |3 < 4 y< 3dJ y /la £ 4pJ> 

-» |3d^«34 7 ,) + hi/ 



(20) 



where O 2p states is integrated out. Although this exci- 
tation process is similaplo that of the Raman scattering 
due to the orbital wave,a the current operators are differ- 
ent between the two. Furthermore, in the Raman scatter- 
ing for the conventional light, the polarization of the light 



A. L-edge case 

The transition probability of RIXS for the orbital wave 
is calculated in the electron-photon coupled system with 
the Hamiltonian, 



H — H e i e + H p h + H e - P 



(21) 



The first term is the electronic one for the Mn 3c? and 
Mn 2p orbitals which includes the following three terms: 



with 



H e le — H 3c [. + H-2p + H 3 d-2p ; 
H 3 d = £d ^ d\ la <jdi~,eo + E d ll<r d jl'c 

+ rii(3d 7eT ) ni(3d lel ) 

ile 

+ U' ^ ri,j(3d 7e+CT ) nj(3d 7e _ CT /) 

+ H e - t + H t -t ; 



(22) 



(23) 



and 



Pi-ycr 



^y(3d 7e ,2p 7 )n i (3d 7e )n w (2p 7 ) . (24) 



*787 



n.i(3c? 7CT ) and nhi{2p ia ) are the number operators for Mn 
3d electron and Mn 2p hole, respectively. V(3d lg , 2p 1 ) is 
the Coulomb interaction between Mn 3d and 2p and is 

given by ^(3d 7e , 2p 7 ) = F (3d, 2p) - 4F 2 (3d, 2p) cos(e + 
The second and third terms in Eq. ( pl| ) de- 
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scribe the photon part and electron-photon interaction, 
respectively, and a sum of them is given by 

Hph + H e -p = ^ cjfe ( a fcA a fcA + 9) ' 

feA 

-E^f ^-e^e-^ + e^al) , (25) 
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where aL is the photon creation operator with momen- 
tum k and polarization A. In the L-edge case, the current 
operator describes the 2p — > 3d transition as 

0k)a = £ A^J lle<yViaa ^ + H.c. . (26) 

The coupling constant depends on the polarization and 
the orbital as follows, 

42*=4>co S (^+n Q f). (27) 

with (n x , n yi n z ) = (1, 2, 3). A$ is a constant part defined 
by 

e / 3y/Z\ I -6z 2 x 2 x 2 z 2 x*\ , no . 

* = ^(^rj(— + TrVr-^v P ) (28) 

where (A) = / drR 2p (r)AR 3d (r). 

Being based on the Hamiltonian, the scattering matrix 
for the orbital wave is derived. The initial and final states 
are described by H^d (Eq.p)) where the orbital wave is 
the eigen-mode. The intermediate state is described by 
the single site term in H e i e (Eq.(p^)) and the inter-site 
effects are neglected. The explicit form of the inelastic 
part of the scattering matrix is obtained as 

(S)f i = -2m8(Ef-E i ) 

X -p- V V(efc 1 A 1 )q(efc 2 A 2 )q 

/cifc 2 AiA 2 v a i 

x e i(^ 1 -o; j!2 )t-»(fe* 1 -fe* 2 ).fl i7 . ajD _ ( 2Q ) 

-E^r/) is the energy at the initial (final) state in the 
electron-photon system. In Appendix A, we discuss the 
elastic part which is a probe to detect the orbital order- 
ing. In the above formula, the orbital excitation by x-ray 
is represented by the operator, 

ha = {A ia l A\ ba Ti- + A ib l A\ aoi Ti + ) , (30) 

where subscript a(b) indicates the occupied (unoccupied) 
3d orbital at site i. D is the energy denominator in the 
perturbation given by 

D = -lE®- 1 + ±EM- 1 . (31) 

Two terms correspond to the different intermediate states 
where two electrons occupy the different e g orbital with 
the triplet and singlet spin states, respectively (see Ap- 
pendix A). It is mentioned that when the wave func- 
tion is chosen to be real, in Eq.(|30|) is expressed by 
T x , and it is in contrast to the neutron scattering where 
both S x and S y appear in the scattering matrix. This re- 
markable feature brings about the selection rule between 
the reciprocal lattice vector and the mode of the orbital 
wave as we will show later. The transition probability 



is obtained by the conventional golden rule. By rewrit- 
ing the pseudo-spin operator in Eq. (|3Fj) by the Holstein- 
Primakoff bosons, we obtain the transition probability in 
the L-edge case: 

W = CY,\F^(q,G)\ 2 

Gqp, 

x (% 1 _g 2+f+ c3^ + ^ ) )(l + ^ ) ) 

+ h 1 -- k2 - ?+ G^~^ ) H^ ) ) , (32) 

with C = (2n/fl y /uJk^uJk^) 2 . is the number of the 

orbital wave of mode /i and G is the reciprocal lattice 
vector. Since the one orbital wave process is possible in 
the present case, the transition probability shows a sim- 
ilar form to that for the one magnon neutron scattering. 
Difference in the present case from the magnon scattering 
reflects on the generalized structure factor: 

(33) 

where V u ^(q) and W vfl {q) are the coefficients in the Bo- 
golyubov transformation which connect the operator for 
the ^-th ion to that for the /i-th eigen mode as 

a vq = ^cf?)^ + w ^^(q) a t- q ■ ( 34 ) 

mu 

It is noticeable that in the above formula, G-, q-, and a- 
dependence of the structure factor is mainly dominated 

by the factors e lG ' r " , V vfJi (q) + W vll (q) and Alfy a , respec- 
tively. 

In Fig. 5, numerical results of the structure factor for 
the spin-F case is presented. The orbital state is chosen 
as (6 A /6 B ) = (§/§ +7r). G = (h,k,l) for the fee lattice is 
chosen as h+k+l = even (Fig. 5(a)) and odd (Fig. 5(b)), 
respectively. It is clear that the (+)- and (— )-modes are 
observed separately in the odd and even cases, and their 
weights are proportional to (cosh 9q + ^ + sinb.9 q + ^) 2 and 
(cosh#g + sinh 9 q ^) 2 , respectively. This selection rule 
is originated from the following facts: 1) when the wave 
functions of the e g orbitals are chosen to be real, the scat- 
tering matrix in Eq.(p9|) and ( f30| ) is represented by T x , 2) 

the absolute value of the coupling constants |A„aaA^ | 2 
in the two sublattices is the same in the orbital order 
(6 a 1 6 A +tt)) 3) the relative motion of the pseudo-spin 
in (+)- and (— )-modes is opposite, and 4) there exists 

the factor e 1 » . It is shown that this selection rule does 
not depend on the phase of the wave functions of the 
e g orbitals. It is highly in contrast to the case of one 
magnon neutron scattering where the correlation func- 
tion of S y breaks the selection rule and the (+)- and 
(— )-modes are also observed in the h + k + l = even 
and odd cases with weights of (cosh^ 4 ^ — sinh0 q + ^) 2 
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and (cosh^g — sinh6>g ') 2 , respectively. The intensity 
of the structure factor monotonically increases with de- 
creasing the energy and it diverges at the T-point. This 
^-dependence comes from V vfl (q) + W ufl (q) which is pro- 
portional to (coshf?^ 1 ' + sinnf?^) 2 . The divergence at 
the T-point is originated from the divergence of the x 
component of the staggered orbital order parameter due 
to the rotational symmetry in the (T z , T x )-plane. In the 
figure, the intensity in the z-polarization is larger than 
that in the x-polarization and the ratio does not de- 
pend on q. It is interpreted that the orbital order is 
close to (3d 3x 2_ r 2 /3e? a 2_ z 2) where A a_ z a x is zero due 
to the symmetry. In Fig. 6, the structure factor in 
the spin-A case with the x-polarization is shown. The 
four modes of the orbital wave, that is, (+±)- and (— ±)- 
modes are mainly observed for h + k=odd, Z=odd (even), 
and h + fc=even, Z=odd (even), respectively. Strictly 
speaking, this selection rule is satisfied in the case of 
J2 = and in general, (++)- and (- -)-modes and (+- 
)- and (-+)-modes are mixed each other. The intensity 
does not diverge at the T-point, unlike that in the spin-F 
case, since the rotational symmetry in the orbital space is 
broken by the anisotropic spin structure. The intensity is 
larger in the presented z-polarization case in comparison 
with that in the x-polarization one, although the results 
in the cc-polarization are only shown in the figure. It is 
attributed to the type of the orbital order which is close 
to (3d 3x 2_ r 2 /3d y 2_ z 2) rather than (3d 3z 2_ r 2 /3d x 2_ y 2). 

B. K-edge case 

In the K-edge case, we derive the effective Hamiltonian 
for the electron-photon system and calculate the transi- 
tion probability by using the Hamiltonian. The electronic 
part of the Hamiltonian, where Mn 3d, Ap and Is orbitals 
are introduced, is given byE3 

H e le = Hq + H 3 d + H 3 d-4p + His-sd + Hl s -ip , (35) 

with 

H a =S P J2 PtyvPvyc +£sJ2 ' ( 36 ) 

if a ia 

where Si a and Pi ia are the annihilation operator of Mn Is 
and Ap electron with spin a and orbital 7, respectively. 
The last three terms in Eq.(|3|) describe the Coulomb 
interactions: 

H 3 d-4p + Hi s -3d + Hi s -ip 

= ^ V{id le ,Ap 1 )n i {id le )n i {Ap 1 ) 

if 91 

+ E V(3d, ls)n,i(3d le )n hl (ls) 

ife 

+ V(4p, ls)m(4py)n hi (ls) . (37) 

if 



In the present case, the current operator in H e - p describes 
the Is — > Ap transition asll3 

(Jfc)a - E ^ K) PL^ m + H.c . (38) 

ia 

A^ K ' is the coupling constant defined by 

A {K) = df^JiOHVcOMO . ( 39 ) 

which is derived by the dipole approximation and does 
not depend on the polarization a. 

The scattering process in the K-edge case includes the 
electron transfer as shown in Fig. 4(b) and is derived 
by the fourth order process with respect to H e - P and H t 
(the transfer term in H 3 d (Eq.(|23])). Therefore, in order 
to avoid complications due to the higher order perturba- 
tion, we derive the effective Hamiltonian by the following 
procedure. In the initial and final states of the scatter- 
ing, the energy of the doubly occupied e g state is higher 
than that of the singly occupied one by the order of the 
Coulomb interaction between 3d electrons. In the inter- 
mediate state, on the other hand, due to the screening 
effects, the energy of the singly occupied state is higher 
that that of the doubly occupied one by the order of (the 
screening potential)- (the Coulomb interactions between 
3d electrons and 3d and Ap electrons) . By excluding these 
two states which have the higher energy eigen values, we 
derive the effective Hamiltonian by the second order per- 
turbational calculation of Ht and/or H e - P as follows, 

H = He-e + H core + H e -p + Hph- (40) 

The detailed derivation is presented in Appendix B. The 
first term was introduced in Eq.(^|) and the second term 
describes the electronic states where the Is core hole and 
the Ap electron always exist. The third term is recognized 
as the effective electron-photon interaction given by 

i8 mkXoi / yy' n f" aa' " 

x (Pmdly^niiSd^di+sy'a') 

x jiaka k x(e k x)ae~ iu '' t + H.c. . (41) 

It implies the electron excitation from Is to Ap accompa- 
nied with the electron transfer from i + 5 site to i site, 
where the operator P m projects out the m-th doubly oc- 
cupied state. is the energy de nom inator derived in 
the perturbational calculation (Eq.(|Bl|)) and (jik)a is the 
i-th component in Eq.(j38|). 

The scattering matrix is calculated by the perturbation 
of the effective interaction H e - P and is given by 

(S) fi = -2mS(E f - E t ) (42) 

X— V y^(efc 1 A 1 )a(efc 2 A 2 )q 

" l U X\ *JUk^k2 

/cife2AiA 2 Ct v «7 



G 



x e 



i(uj kl -a>fe 2 )i— i(fci— k 2 )-Ri 



(43) 



^i m) and are denned in Appendix B. J^j is the 
main part of the formula and is expressed by the spin 
and orbital pseudo-spin operators in i and its nearest 
neighbor i + S sites. It is noticeable that is de- 

scribed by T x and T z (see Appendix B), as well as that 
in the L-edge case in Eq.([30|). It describes changes of 
the spin and orbital states in these sites. The scattering 
processes, where the one and two orbital waves are con- 
cerned, are included in 1^ and these processes occur in 
the same order. This is in contrast to the conventional 
magnon Raman scattering where the two magnon sector 

is dominant. This is attributed the fact that £Q 7^ (5 77 <. 
Therefore, we consider the processes for the one orbital 
wave excitation in schematically shown in Fig. 7. 
In this approximation, the transition probability is ob- 
tained as the same form as Eq. (J3^) with the structure 
factor replaced by 

vm 

+ e i ^-^{V vl ^q) + W v . tl {q))B < ^) . (44) 

The index v' indicates one of the nearest neighboring Mn 
ions of the i^-th ion and B^) and Cffl are represented by 
the spin correlation function defined in Appendix B. The 
two terms in the parenthesis in Eq. ( f44| ) corresponds to 
the orbital excitation in v site and its nearest neighboring 
v' site, respectively. As a result, the factor e t ^'^ Vl '~ r "'^ 
gives the additional (^-dependence in the structure factor. 

In Fig. 8, the numerical results of the structure factor 
in the spin-F case are presented. The parameter values 
are chosen as u>o/io = 0.5, W/t n = —6, U/t n — 5, J/t = 
1.0, and T/t — 0.5 (see Appendix B), although the qual- 
itative results do not depend on the parameter values. 
h + k + Z=odd and even cases mainly correspond to (+)- 
and (— )-modes, respectively, although the both intensi- 
ties are strongly mixed each other. The q-dependence of 
the intensity is not monotonic, unlike that in the L-edge 
case shown in Fig. 5. Especially, around the T-point, 
the intensity gradually decreases and becomes zero at 
the point. This is attributed to the following unique 
excitation process in the K-edge case. With using the 

relations = and E,<f w *W = °> in the 

uniform spin structure, C vv < = and ~Y^ V , B vv i = 
are derived. These relations suppress the divergence of 
the factor V^q) + W^cf) at the T-point in Eq. (44), so 
that the intensity becomes zero at the point. It implies 
that at the T-point the several orbital excitation pro- 
cesses by x-ray are canceled out each other. The other 



characteristic cf-dependences of the intensity, such as, no 
intensity along the T — L direction, are interpreted by the 
factor e i y'( r "~ r v'\ The polarization dependence of the in- 
tensity comes from the Coulomb interaction V(id 7 ,4p a ) 

included in r 7 ™\ although it is not remarkable in com- 
parison with that in the L-edge case (Fig. 9). Around the 
X-point, the slightly larger intensity is shown in the x- 
polarized case in comparison with that in the ^-polarized 
one and it is interpreted as follows. In the present orbital 
order, a dominant excitation process around X-point is 
that in Fig. 5(b), where the orbital states in i and j 
sites are 9 = | + it and f , respectively. Since the orbital 
at the i-site, that is, ^(3d3 Z 2„ r 2 + 3d x 2 -3,2), shrinks in 
x-direction, the Is — > Ap x excitation shows a larger in- 
tensity than that of Is — > Ap y ( z ) excitation due to the 
Coulomb interaction. Therefore, the intensity is slightly 
larger in the z-polarized case. In Fig. 10, the results 
in the spin-A case with the s-polarization are presented. 
The relation between the reciprocal lattice vector and the 
four excitation modes is the same as that in the L-edge 
case. However, the two modes are mixed and the mix- 
ing becomes remarkable around the M -point due to the 
factor e l Q'( rv ~ r »'\ The intensity of the structure factor 
does not show a strong polarization dependence of x-ray 
in comparison with that in the L-edge case, since the po- 
larization dependence is originated from that in as 
explained above. 



IV. SUMMARY AND DISCUSSION 

In this paper, we study the orbital excitation in the 
orbital ordered manganites and the method to observe 
the excitation by using the resonant inelastic x-ray scat- 
tering. At first, we emphasize the noticeable correlation 
between the static spin structure and the orbital exci- 
tation. The anisotropic spin ordering brings about the 
anisotropy in the orbital space and causes the gap in the 
orbital excitation. In this case, the interaction between 
spin wave and phonon seems to be weak, since the gap 
in the orbital excitation is of the order of t^/U . This 
situation is expected to be realized in LaMnOs where 
the orbital ordered temperature accompanied with the 
structural phase transition occurs To — 780A' and the 
A-type antiferromagnetism is realized at Tm = 140A.H3 
Although the coupling with the lattice contributes to the 
gap, it is predicted that the nature of the orbital ex- 
citation is changed at Tjv and below the temperature 
the gap becomes more remarkable. On the other hand, 
in the spin-F case, the orbital excitation becomes gap- 
less and the low lying excitation affects the thermody- 
namic and transport properties. In this case, the disper- 
sion relation is modified by the interaction between spin 
wave and phonon. The specific heat which is propor- 
tional to T 3 is derived by the linear dispersion of the or- 
bital wave, which is distinguished from the contribution 
from the ferromagnetic spin wave. The possible candi- 
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date, where the gapless excitation is anticipated, is the 
ferromganetic insulating phase realized in the low hole 
doped manganites, where the superexchange interaction 
between e g orbitals dominates the ferrom agn etic interac- 
tion rather than the double exchange one.EJO Especially, 
in Lao.88Sro.i2Mn03, the orbital ordering is confirmed by 
the resonant elastic x-ray scattering below Tl = 140K, 
and the static Jahp»Teller distortion is released in this 
temperature rangeJlj Therefore, the low lying excitation 
is expected to be observed without disturbance of the 
lattice distortion. In the spin canted phase among the 
above two concentration region, the anisotropy of the 
spin correlation function decreases and the gap in the 
orbital excitation gradually decreases with doping holes. 
Such kind of control of the orbital excitation is possible 
by changing the spin structure in the magnetic field. 

As a probe to detect the orbital excitation, we propose 
the two possible methods by using the inelastic resonant 
x-ray scattering. In the present case, the dispersion re- 
lation is detectable due to the smaller wave length of 
the x-ray in comparison with the conventional light (the 
wave lengths in the L- and K-edge cases are about 2.1A 
and 19. 6A, respectively). In the L-edge case, the orbital 
excitation is directly caused by 2p — > 3d transition and 
the single orbital wave without the spin flipping is only 
brought about. In the spin-F case with orbital (0/6 + tt), 
each mode of the orbital wave is separately observed by 
choosing the reciprocal lattice vector G. This selection 
rule between the reciprocal lattice vector and the modes 
of the orbital wave is mainly due to the fact that the 
inelastic part of the scattering matrix is represented by 
the x component of the pseudo spin operator. Since the 
large scattering intensity is expected in the low energy 
region near the orbital superlattice reflection point, it 
may be observed without disturbance of the large fun- 
damental elastic peak. The polarization dependence is 
remarkable, since it directly comes from the coupling 
constant between electron and photon. In the K-edge 
case, where the elastic scattering has already been used 
to observe the orbital ordering, the orbital excitation oc- 
curs through the screening process of the core hole. This 
process through the charge transfer from one of the near- 
est neighboring Mn ions is more favorable in manganites 
rather than in cuprates and nickelates, since the man- 
ganites are closer to the boundary between the charge 
transfer insulator and the Mott one. The one and two 
orbital wave excitations are possible in the same order, 
unlike the conventional two magnon Raman scattering. 
The orbital excitation occurs in the site where the x-ray 
excites the Is electron and its one of the nearest neigh- 
boring Mn sites. Owing to the unique excitation process, 
the structure factor shows the strong momentum depen- 
dence and in particular the intensity vanishes around the 
T-point. Through the measurement of the characteris- 
tic momentum dependence, it is possible to identify the 
orbital sector in the inelastic component and the mecha- 
nism of the orbital excitation proposed in this paper. 
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APPENDIX A: ELASTIC SCATTERING IN THE 
L-EDGE CASE 

In this appendix, we introduce the elastic x-ray scat- 
tering in the L-edge case as a method to detect the or- 
bital order. We focus on the polarization and orbital 
dependence of the atomic scattering factor in the orbital 
ordered phase. Once the dependence is revealed, the scat- 
tering intensity in the several orbital ordering is obtained 
by the formulas in Ref. |l4|. The atomic scattering factor 
is calculated by the perturbational calculation presented 
in Sect. Ill A and given by 

fiaa = J a f EM-' 

+ i4^J 2 ^ (s) " 1 + i4^i 2 ^ (d) " 1 - ( A1 ) 

As well as in the inelastic case, we approximately set 
up the state and the energy for the intermediate states 
in the scattering. We assume that the states consist of 
the following three states: \t2p a ), \s2p a ), and |d 7 2p Q ), 
where |s)(|t)) represents the state where two e g elec- 
trons occupy the different orbital with the spin singlet 
(triplet) and |eL) is the doubly occupied state in the 7 
orbital. These energies are approximately estimated as 



= U'-J+A+2F (3d,2p)-Lu k +iT with A = e d -e 
J5« = U' + J + 2 J H + A + 2F (3d, 2p) - uj k + iT and 
E^ = U + 2 J H + A + 2V(3d n , 2p a ) - uj k + iT, where 
r is the damping constant. Because E^ is the smallest, 
\t2p a ) dominates the L-edge. It is noticed that the en- 
ergy is independent of the occupied 3d orbital and 
the polarization of x-ray, because two e g orbitals are oc- 
cupied in the intermediate state and three 2p orbitals 
are occupied in the initial and final states. On the other 
hand, its intensity in Eq.(Al) is proportional to |j4 7 ^_a| 2 
which clearly depends on the occupied orbital and the po- 
larization as shown in Eq.(p7|). It is concluded that the 
anisotropy of the scattering factor is originated from the 
coupling constant between electron and photon, unlike 
that in the K-edge case where the Coulomb interaction 
between 3d and 4p gives rise to the anisotropics For ex- 
ample, In the antiferro-type orbital ordering, where two 
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kinds of the orbital sublattice exist, the scattering inten- 
sity at the orbital superlattice reflection is proportional 
to |(|.A 7A _ Q | 2 — |^4 7b _q| 2 )| 2 , where 7.4(B)- is the unoccu- 
pied orbital state in the A(B)-orbital sublattice. 



APPENDIX B: DETAILED DERIVATION AND 
FORMULAS IN K-EDGE CASE 

In this appendix, we present the detailed derivation 
and formulation of the effective Hamiltonian and transi- 
tion probability in the K-edge case which are introduced 
in Sect. Ill B. 

In derivation of the effective Hamiltonian, the unper- 
turbed states are classified by the number of electron 
of Mn 3c? and 4p orbitals at each site as \n(3d), n(Ap)) 
and these states and energies are approximately set up 
as follows. The energy in |1,0) is defined as the origin 
of the energy. |1,1) consists of |3rf^,4pLls) and its en- 
ergy is dinned as E n = A + V{3d 1 ,Ap a ) + V(3d, Is) + 
V(4p, Is) — LOk with A = ep — e s . |2,0) includes the 
three states represented by \t), \s) and \d), as we ex- 
plained in Appendix A, with the energies E_\q = U' — J, 
E$ = U' + J + 2J H and E$ = U + 2J H , respec- 
tively. In the same way, |2, 1) consists of the three 
states: \t 4p*ls), \s 4p* Is) and |d 7 ip^ls) with en- 
ergies E$ = U' - J + A + 2F (3d,4p) + V core - cj fe , 
E ( 2 f = U' + J + 2J H + A + 2F (3d,4p) + V core - iu k 

and E$ = U + 2J H + A + 2V{3d 1 Ap a ) + Vcore - w*, 
respectively, with V core = 2V(3d, Is) + V(ip, Is). In 
all states, the ti g spin is assumed to be parallel to the 
e g one and the spin dependence of the core hole in- 
teraction are neglected. Among these parameters, we 
choose the independent parameters wq = AF% (4p, 3d) , 
W = U + 2 J + V(ls, 3d) + F (4p, 3d) - 5 J and U = 
U-3J with U' -J + A + 2F (3d, 4p) + V core - u k = 0, 
V(ls,3d) = V(ls,4p) and J = J H . 

The matrix element of the effective Hamiltonian (H)i n 
is derived by the conventional canonical transformation. 
The perturbational processes are characterized by the set 
of the initial, intermediate and final states: |m), \n)). 
In the calculation, \n(3d), n(4p)) = {|1,0), |2, 1)} are in- 
cluded explicitly and {|2, 0), |1, 1)} are introduced in the 
virtual senses. H 3< i in Eq.(|40|) corresponds to the pro- 
cess (\l),\m),\n)) = (|1,0),J2,0),|1,0». The effective 
electron-photon interaction H e - p in Eq. (041) is derived by 
(|2,1),|2,0)|1,0)) and (|2, 1), |1, 1), |1, 0))7ln this term, 

1 \oi is given by 

r^ = 1 -((E^-E 11 y 1 + (E 10 -E 11 y 1 

- E^y 1 + (e w - E^y 1 ) • (bi) 

The second term H core in Eq.(|40|) is derived by the pro- 
cess (|2, 1), |1, 1), |2, 1)) and its 0-th order term of the 



perturbation is used in the calculation of the transi- 
tion probability. The processes (|1, 0), |1, 1), |1, 0)) and 
(|2, 1), |2, 0), |2, 1)) imply the elastic off-resonance scat- 
tering and the higher order x-ray scattering, respectively, 
and these are irrelevant in the present scattering. 

In the scattering matrix in Eq. ( |43|) , is the energy 

denominator originated from the second order perturba- 
tional process of H e - p and it is represented by 

D^ = {E^ ) -E 10 + iTy 1 . (B2) 

in Eq.(p3|) describes the changes of the spin and 
orbital states defined as 

4jl = (4 + Si ■ &j) 

\ ba 'ij T 1 act 1 ij ' 1 aa L ba ' ij ) ' \ D ^J 

and 

$ = 2(i-3-£,)(r«V+ r i?#). (B5) 

with 

4 = "a (t?fn ja + t% 2 n jb + 2t%t$T jx ) , (B6) 
= n la {t\fn ja + iff n i6 + 2t%t%T jx ) , (B7) 

C c\ f ,aaj.barp , ,bbj_abrp , tyi.aaj.bbrp rp \ 

I ij ^^ij L ij 1 ix fl ja ~r L ij L ij 1 ix fl jb ~r ^t-jj l i j±i x ±j x } , 

(B8) 

r* = n ia {t?fn ja + t$ 2 n jb + 2t%t$T jx ) , (B9) 

t? = n lb (t\fn ja + t™n jb + 21%1%T jx ) . (B10) 

The coefficients B { ^ ] and in Eq. @ are ob- 

tained by taking the average for the spin part as 

c { 3 = -(- + $,- M)rMc<C' . (Bii) 
c% - (\ 0u ■ i))rWrWcft > (B12) 

C^J = , (B13) 
S$ = (| + & ■ S^y^HlUt, , (B14) 







and 



- <& • &>)r$ a C'*& ■ (Bi6) 
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Figure captions 

FIG. 1. The mean field phase diagram at zero temper- 
ature. F, A, C, and G imply the ferromagnetic struc- 
ture, and the layer-type, rod-type and NaCl-type anti- 
ferromagnetic structures, respectively. 

FIG. 2. The dispersion relation of the orbital wave in 
the spin-F case. The Brillouin zone for the face-centered 
cubic lattice is adopted. The orbital states for the two 
orbital sublattices are denoted by (6 a/ 6 a + tt)- 

FIG. 3. The dispersion relation of the orbital wave in the 
spin-A case. The Brillouin zone for the tetragonal lattice 
is adopted. J2/J1 = 0.5 which corresponds to the orbital 
state (6a/ — 6a) with 6a = 1-34. \i denotes the mode of 
the orbital wave. 

FIG. 4. The schematic picture of the orbital excitation 
processes in (a) the L-edge case and (b) the K-edge case. 
j denotes one of the nearest neighbor sites of i. The 
broken arrows indicate the incident and scattered x-rays. 

FIG. 5. The square of the structure factor for the 
orbital wave with spin-F in the L-edge case: II — 
\F^ a (q,G)\ 2 /(D\A \ 2 ) 2 . (a) h + k + l=even, and (b) odd 
where (h, k, I) is defined in a unit cell in the face-centered 
cubic lattice. The orbital state is chosen as (§/f 
The shaded circles and squares mean the x and z polar- 
ization cases, respectively. 

FIG. 6. The square of the structure factor for the 
orbital wave with spin-A in the L-edge case: II = 
\F^ a (q,G)\ 2 /(D\A Q \ 2 ) 2 . (a) h + fc=cven and l=odd, and 
(b) h + fc=odd and Z=odd, where (h, k, I) is defined in a 
unit cell in the tetragonal lattice. J2/J1 = 0.5 which cor- 
responds to the orbital state (6a/ — 6a) with 6a = 1-34. 
The polarization is parallel to the rr-direction. The open 
and filled circles mean (a) the (+— )- and ( — h)-modes, 
and (b) the ( )- and (++)-modes, respectively. 

FIG. 7. The schematic picture of the excitation process 
of the one orbital wave in the K-edge case, j denotes 
one of the nearest neighbor sites of i where Is — > 4p 
transition is caused by x-ray. The two levels indicate the 
degenerate e g orbitals. The straight and broken arrows 
indicate the first and second electronic processes in the 
scattering, respectively. 

FIG. 8. The square of the structure factor for the 
orbital wave with spin-F in the K-edge case: Ik = 
\F^ a (q,G)\ 2 /(2t^ 1 \A( K 1\ 2 ) 2 . (a) h + k + l=even, and 
(b) odd where (h, k, I) is defined in a unit cell in the 
face-centered cubic lattice. The orbital state is chosen as 
(§/§ + 71 ")- The polarization is parallel to the x-direction. 
The open and filled circles mean the (— )- and (+)-modes, 
respectively. 

Fig. 9 The polarization dependence of the square of the 
structure factor for the orbital wave ((+)-mode) with 
spin-F in the K-edge case, (a) h+k+l=even, and (b) odd 



where (h, k, I) is defined in a unit cell in the face-centered 
cubic lattice. 

FIG. 10. The square of the structure factor for the 
orbital wave with spin-A in the K-edge case: Ik — 
\F l , a (q,G)\ 2 /(2to 1 \A^\ 2 ) 2 . (a) h + fc=evcn and Z=odd, 
and (b) h + fc=odd and Z=odd, where (h, k, I) is defined 
in a unit cell in the tetragonal lattice. J2/J1 = 0.5 
which corresponds to the orbital state (6a/ — 6a) with 
6 a = 1.34. The polarization is parallel to the x-direction. 
The open and filled circles mean (a) the (H — )- and ( — (-)- 
modes, and (b) the ( )- and (++)-modes, respectively. 
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